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Abstract. We deduce an overcomplete free energy functional for D = 1 par- 
ticle systems with next neighbor interactions, where the set of redundant vari- 
ables are the local block densities Qi of i interacting particles. The idea is to 
analyze the decomposition of a given pure system into blocks of i interacting 
particles by means of a mapping onto a hard rod mixture. This mapping uses 
the local activity of component i of the mixture to control the local association 
of i particles of the pure system. Thus it identifies the local particle density of 
component i of the mixture with the local block density Qi of the given system. 
Consequently, our overcomplete free energy functional takes on the hard rod 
mixture form with the set of block densities Qi representing the sequence of 
partition functions of the local aggregates of particle numbers i. The system 
of equations for the local particle density q of the original system is closed via 
a subsidiary condition for the block densities in terms of q. Analoguous to the 
uniform isothermal-isobaric technique, all our results are expressible in terms 
of effective pressures. We illustrate the theory with two standard examples, the 
adhesive interaction and the square-well potential. For the uniform case, our 
proof of such an overcomplete format is based on the exponential boundedness 
of the number of partitions of a positive integer (Hardy- Ramanuj an formula) 
and on Varadhan's theorem on the asymptotics of a class of integrals. We also 
discuss the applicability of our strategy in higher dimensional space, as well 
as models suggested thereof. 



1. Introduction 

The era of the class of exactly solvable D — 1 particle systems with next neighbor 
interactions in classical density functional theory (for an overview see, e.g., JlO[) 
began in 1976, when Percus derived the free energy density functional for hard rods 
in an arbitrary external potential held j8|. Soon afterwards, he presented the inverse 
operator format B . This formalism not only allowed him to simplify the derivation 
of the hard rod functional, but also led him to the solution of the sticky core 
model, in which the interaction has a 0-range attractive component. Moreover, he 
showed that — analoguous to the uniform case — the adhesive interaction is the only 
potential with an attractive component, where an explicit solution is possible. The 
generalization of the pure hard core result to an additive mixture of hard rods was 
established by Vanderlick, Davis, and Percus in 1989 jL3|. Already a few years later, 
Brannock and Percus constructed along Wcrthcim's theory of local association (for a 
stringent review of Wertheim's approach see Jl(i[| ) a free energy density functional for 
multicomponent systems with arbitrary next neighbor interactions [|| . One of the 
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specific examples treated was the one-component case of the prototypical square- 
well interaction, called Herzfeld Goeppert-Mayer (HGM) system || . In abstracting 
from this Wertheim local association free energy density functional format, Percus 
invented in 1997 the overcompleteness technique p[ . In this paper, we present a 
particular realization of this strategy. 

The essential idea of an overcomplete description ]ll[ of thermodynamic func- 
tionals is to introduce a set of additional variables (over the basic local particle 
density g) such that with respect to a direct evaluation the overcomplete format 
has a simpler formal structure and is physically more transparent. Thus, apart 
from the purely mathematical problem of deriving such a representation from first- 
principles, the second difficulty [ pT[ is to give a physical interpretation for the 
redundant variables. 

In our approach we use an observation of Kierlik and Rosinberg |7j to overcome 
the latter difficulty. Their idea, although stated originally only for nonuniform 
dimerizing hard rods (for the uniform case see, e.g., |5[), consists in interpreting a 
D = 1 particle system with next neighbor interactions as composed of blocks (clus- 
ters [||, molecules || Q, or superparticles Q]) of i interacting particles (monomers, 
dimers, trimers, ... ) and to use the corresponding local block densities Qi as the 
set of additional variables. Here two particles are understood to interact, if in the 
standard decomposition of the Boltzmann factor into a reference hard core and an 
attractive part, ( |l.3|) , the interaction between the two particles is described by the 
latter one. Since the blocks of local aggregates interact by construction via hard 
core potentials, the free energy functional takes on the simple hard rod mixture 
form. In such an overcomplete description the physical significance of the redun- 
dant variables is given at the outset. 

So our program is to analyze a given pure system with arbitrary next neighbor 
interactions by means of a mapping onto a hard rod mixture. This mapping uses 
the local activity of component i of the mixture to represent the local aggregates 
of i particles of the pure system and so it identifies the local particle density of 
component i of the mixture with the local block density Qi of the original system. 
Consequently, with the free energy of hard rod mixture form, the set of local block 
densities represents the sequence of partition functions of the local aggregates with 
particle numbers i. We close the system of equations for the local particle density 
via a relation between g and the block densities inferred through a trivial identity 
of functional analysis. 

Thus already implicitly formulated by Kierlik and Rosinberg 0, they settled 
this general program only for the special case of dimerizing hard rods. They argued 
that, although this technique can be extended to the more general case of sticky 
cores, where blocks of any number of particles coexist, the possible difficulty of this 
method lies in the definition of blocks for arbitrary next neighbor interactions. So 
they proceeded with Wertheim's theory instead. They (i) showed that Wertheim's 
first-order thermodynamic perturbation theory (TPT1) is exact for dimerizing hard 
rods as well as for sticky cores and (ii) derived density functionals for these systems 
with the simple hard rod mixture form through TPT1. Later on, the exactness of 
the TPT1 also for sticky core mixtures was established by Brannock and Percus 

Actually, already Brannock and Percus used the block interpretation of Kierlik 
and Rosinberg |Q to construct an overcomplete free energy functional. But only 
formally, as a purely technical intermediate step in their proof, because their aim 
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was to infer density functionals that take on the Wertheim local association form. 
However, in such an overcomplete description the remarkable hard rod mixture form 
of the free energy is obscured. As already stated, again Kierlik and Rosinberg 
first revealed that at least for dimerizing hard rods and for sticky cores, functionals 
of the simple hard rod mixture form underly the Wertheim functionals. 

Similarly, also with respect to the relative density format of Percus , another 
earlier and formally alike equivalent overcomplete description, where the redundant 
variables are partition function densities, our realization has a simpler formal char- 
acterization (in terms of a functional of hard rod mixture form) and is physically 
more transparent (with the block densities as the basic elements of the theory). 
Moreover, all our main results possess higher dimensional analogues, so that our 
overcomplete description immediately suggests D > 1 models as well. 

Formally, we start from the (configurational) isothermal measure (or the Gibbs 
measure thereof) for a particle system on the line R with next neighbor interactions, 
given by the density 

n— 1 n 

\ [ e{q k+1 ,q k )l[exp[-pU(q k )} (1.1) 

k=l k=l 

with respect to the usual Lebesgue measure. Here, n denotes the number of parti- 
cles, their positions being q k el, U is the external potential, and the translation- 
invariant interaction V is represented by the Boltzmann factor (at reciprocal tem- 
perature /3) 

e(y, X ) = h p[ - my - x)] i0Ty ^ X > (1.2) 
I for y < x. 

We assume that V has a hard core of diameter a > 0. Hence if we restrict ourselves 
to the next neighbor case, the range of V is 2a, i.e., V(x) = 0, x > 2a. Furthermore, 
we introduce the standard decomposition (for a related analysis based on a different 
decomposition see, e.g., Q) 

e(y,x) = h a (y-x) + f(y,x) (1.3) 

with h a the Heaviside function shifted to the right by a and f(jj,x) = h a (y — 
x)[e(y, x) — 1], known as the attractive part of the Boltzmann factor. 

First, we will treat for the sake of completeness the uniform case within our 
overcomplete description. This will also allow us to introduce the concept most 
clearly. Our account will be inductive. We are going to develop the theory along 
two standard examples, the sticky core model and the HGM system. Second, we 
will lift our overcomplete format to the general nonuniform case and illustrate the 
results with the two examples already examined for the uniform case. Finally, we 
will indicate how the idea extends to higher dimensional space and suggest D > 1 
models thereof. 



2. Uniform case 

In this section we introduce the concept of overcomplete description for the free 
energy of uniform systems with arbitrary next neighbor interactions. We derive 
such a representation from first-principles as n — > oo based on the exponential 
boundedness of the number of partitions of a positive integer n (Hardy- Ramanuj an 
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formula, see, e.g., |ffl]) and on Varadhan's theorem on the asymptotic behavior of a 
class of measures (see, e.g., 

Formally, we assume U(x) = 0, x G M. Then the total isothermal measure 
on (0, l) n , under the condition that the particles be reflected elastically at the 
endpoints of the interval (0,1) (the volume |(0, Z)| = I is called the length of the 
system), becomes 

~ n — l 

exp[-n/3f n (l)]= T[e(q k+1 ,q k )d(qi,...,q n ), (2.1) 

where we have introduced /„ known as the (configurational) n-particle free energy 
per particle. 

For the sake of explicitness and to make contact with a previous paper [Q , we 
start our presentation with the sticky core model. 

2.1. Adhesive interaction. Adhesive particle systems are systems composed of 
hard rods with an additional attractive interaction on the boundary of the particles, 
which allows for local association. Thus an adhesive system may be represented as 
a limiting case of a HGM system. In such systems, the interaction is 

^)=( + ;° (2.2) 

{- £ X[0A)( 2L r) for x G [a,oo), 

where d is the width of the well (note that we have restricted ourselves to d < a), 
£ the interaction strength, and xa denotes the characteristic function of the set A. 



Then — on inserting ( |2.2| ) into ( [L5| ) via (1.2) — adhesive systems are characterized 
by the limit 

lim f(y,x)\ de0£ ^ x = XS a (y- x) (2.3) 

o — *oo 

with S a = h' a , A = "fa, and 7 G M+ being the stickiness parameter of the interaction. 



Hence, due to the O-range attractive component (2.3), our overcomplete descrip- 
tion for the n-particle sticky core model translates to the decomposition into blocks 
(monomers, dimers, trimers, . . . , n-mers) of pure hard core particles with diame- 
ters (a±, . . . , a n ), a,i — ia, and to the introduction of the n block densities as the 
set of redundant variables. Therefore /„ is defined on an enlarged space with n 
additional dimensions. Consequently, the free energy is calculated by means of a 
minimum principle, ( [2.14 ). 



Formally, upon substituting ( J2.3D into (2.1) via ([1.3]), /„ can be interpreted 
as a sum over a binary tree with 2" _1 branches. Consider the contribution of an 
arbitrary branch, i.e., an arbitrary convolution sequence of h a and A<5 a , to this sum. 
For y, x <E (0,1) denote by a block (or cluster Q, molecule ^ |7j, or superparticle 
|2|) bi of i G {1, . . . , n} interacting particles the convolution sequence 

h/ ^ f h a(y~x) fori=l, 

bi{y,x) = < (2.4) 
I Oj_i * A d a (y — x) tor 1 < 1 < n 

with a * b(y, x) = J Q a(y — z)b(z — x) dz. Let fcj be the number of times bi (modulo 
h a at the right boundary) appears in such a branch. For each n > 2 define the set 

An = i k = (fci, ...,k„): ki G {0,...,n} A^ih = n \. (2.5) 
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Fix k € A n . Then, on using the commutativity of the convolution (for another 
application of this elementary property see, e.g., H), we conclude that each of the 



C(n,k) = 



(2.6) 



branches characterized by the vector k contributes the same value to /„. Again 
through the commutativity of the convolution, we can rearrange the blocks accord- 
ing to their length, starting with the group of k\ particles of length a next to the 
origin, so that the total isothermal measure reduces to 



exp[-n/?/„(Z)] = C ( n ' k ) ^' =l(i_1) ' i exp[-n/?/ n , fc (Z)] 

keA n 



(2.7) 



with 



£{exp[-n/3/„ !fc (/)]}(s) = -cxp(as)s 



exp(— ias) 



k, 



l expj- {YJLii-h - I) as] 1 



(2.8) 



where £{exp[— nPf n ^k(l)]}(s) denotes the Laplace transform of exp[— n(3fn,k(l)] 
with respect to I and we used the relative positions 



a^o = ?1, x k = q k+ i - qk for 1 < k < n, x n = I - q„ 



(2.9) 



as well as (2^4). On transforming back, we infer 



exp[-n/?/„(Z)] = /i( n -i)o(0 C ( n > fc M S 



(2.10) 



Finally, let us analyze the leading order asymptotic behavior of /„ for the non- 
trivial case I > (n — l)a. Denote by ir(n) the number of partitions of a positive 
integer n. Since \ A n \ = 7r(ra), it follows by the leading term of the asymptotic series 
of the Hardy-Ramanujan formula (see, e.g., jlj) that as n ~ + oo 



\An\ 



2ttV2 dn 



<,X1,, ^V""24 



24 



(2.11) 



Hence the number of terms in the sum ( 2.1 0| ) is exponentially bounded. Thus one 
may determine the asymptotic behavior of f n by the asymptotic behavior of the 
largest summand in (2.10). This is a standard argument of large deviation theory 
(see, e.g., Q). Therefore, if we use a weak form of Stirling's formula, ln(n!) = 
n In n — n + 0(ln n) , we find as n — > oo 
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and so, on introducing the particle density p 1 = l/n with p < a 1 as n -> oo by 
I > (n — l)a, via (2.11) that 



lim (3f n (p) = lim min { 'S~* — 

K »=1 



In - 1 - In A 1 
n 



In p 



A- 
« — 

71 

8 = 1 



(2.13) 



By the properties of the sets {A„}„> 2 , of the set A = {t = (tx, ■ ■ ■ ) : U > A 
Sfci — l}i an( i 01 tne function on the right hand side (RHS) of ( 2.13 ), we 
finally infer 



lim Pf„(p) = Pfoo(p) = min 

n— *oc c£A 



^ oo 

{E^h^r-i-^- 1 -*) 

^ i — l 



(2.14) 



where /oo is the free energy per particle. This proves (4.2)-(4.4) of [0. 

Fix 7,(ie R + , p e (0, a -1 ). The function appearing on the RHS of ( 2.14 ) is 
convex. Hence it attains its minimum at the unique value c — (ci, . . . ) with 

1 



(i + xoy 

Cl 



1 + Ye 



for i = 1, 
for i > 1. 



(2.15) 



Here, while it is possible to write c directly as a function of p, we preferred to 
introduce for later reference the pressure p via the thermodynamic relation 



f3p{p) = p 2 ^-pfoo{p) 



dp 



with 6 = j3p{p). Thus we have by 



i 



1 + A0 



(2.16) 



(2.17) 



that 



(1 + \0) 



(2.18) 



or explicitly 



Therefore /oo simplifies to 

/3/ooGo) 



1 

2A 



'1 + 



4A 



p 1 — a 



ln(l + A0) -6{p- 1 - a) +ln( 



(2.19) 



(2.20) 



which is the known result. Formula (2.15) is equivalent to the findings (4.5) and 
(4-7) of @. 

As will be shown in Section the remarkable hard rod mixture form of /oo as 
given by (2.14), is an instance of a general result that even holds true in the context 
of nonuniformity 
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2.2. Arbitrary next neighbor interactions. So far, we treated particle systems 
with adhesive interactions. Essentially, the full and explicit analysis was possible 
because of the contact nature of the adhesive potential. When arbitrary next 
neighbor interactions are present, only the minimum of the length of a block hi 
of i interacting particles is ia. The purpose of the present subsection is to establish 
our overcompleteness technique for the general uniform case. 

Formally, our proof of such an overcomplete description for uniform systems with 
arbitrary next neighbor interactions as n — > oo is based on the Hardy- Ramanuj an 
formula, ( [2.11 ), and on an application of Varadhan's theorem on the asymptotic 
behavior of a class of measures (see, e.g., H). 



(L 

by 



Consider first the prototypical HGM system. Then, upon substituting (2.2) into 



via (1.2), the total isothermal measure on (0, l) n is completely characterized 



f{y,x) = \ h a , d (y - x) 



(2.21) 



with h a ^ = h a — h a+ d and A = eJ 3£ — 1. Hence, as in the deduction of our over- 
complete description for the adhesive interaction, we start the proof by analyzing 
an arbitrary branch of the binary tree associated with /„. Let us again decompose 
such a branch into ki (modulo h a at the right boundary) blocks bi of i € {1, . . . , n} 
interacting particles, now defined by the convolution sequence [cf. (2.4)] 



h(y,x) 



h a (y - x) for i = 1, 

\bi-i * Xh a ,d(y - x) for 1 < i < n. 

Then, via (U) as well as ( |2.22| ) and on using the commutativity of the convolution, 
we conclude that 

exp[-n(3f n (l)] = C(n,fc)AS?= 1 ( i - 1 )^exp[-n/3/ n , fe (Z)] 



(2.22) 



(2.23) 



with 



£{exp[-n/3/„,fc(/)]}(s) = i exp(as) s TT < exp ( ias ^ 



1 — cxp(— ds) 



exp(— as) 



[1 — exp(— ds)] 



E?=i(i-i)fe« 



k G A n . 
(2.24) 

Next, let as usual / = np~ 1 , and consider the asymptotic behavior of f n for 
the nontrivial case p^ 1 > a. Since \A n \ is exponentially bounded as n — > oo 
by the Hardy-Ramanujan formula, (2.11), the asymptotics of /„ is governed by 
the asympto tic b ehavior of the largest summand in ( 2.23 ). Thus, by the formal 
structure of ( [2.23] ), the analysis of /„ as n ~ + oo reduces to an asymptotic evaluation 
of f n ,k- If we combine such an evaluation of f n _k as n — > oo by means of Varadhan's 
theorem (see, e.g., Q) with ( |2.12| ), we finally infer for j3 € (0,oo) 



Pfoo(p) 



mm 



In 6 



l-^c,)ln[l-exp(^)] 



(2.25) 



s 
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The case (3 = oo has to be treated separately (see, e.g., Q). 

Fix a, p~ x ,£ ,d, (3 £ R+ such that d < a and p^ 1 > a. Then the function 
appearing on the RHS of (HI) is convex. Hence the minimum is attained at the 
unique value c = (ci , . . . ) with 

1 



{l + A[l-cxp(-0d)]}^ 

Cl 



1 



1 



and thus 



A [1 - exp(-6W)] 



1 



for i = 1, 



for i > 1, 



(2.26) 



i=i 



1 + A [1 - exp(-0d)] ' 
Therefore /oo simplifies for (3 S (0,oo) to 

A/ooOO = - ln{l + A [1 - exp(-0e0]} - 6 (p- 1 - a) + In 
which is the known result. In the sticky limit £ — * oo A d 



(2.27) 



(2.28) 

such that de@ £ — > ja, 



expressions ( 2.28 ) respective ( 2.26 ) reduce to ( 2.2C ) respective ( 2.15 ). An approxi- 
mation of ( 2.26 ) was g iven by He rzfeld and Goeppert-Mayer already in 1934 ||. 
Moreover, through (2.16) and ( [2.25 ), we have the interpretation 



6 = 0p(p), (2.29) 

and so, upon using d2~26| ) in j2~25| ), that the pressure for f3 € (0, oo) is uniquely 
determined by 



fip 



exp((3pd) ( 1 



A 



(2.30) 



1 



Hence, p is a differentiable function with range R+. In particular, we have ^ > 0. 

Finally, since Varadhan's theorem allows us to analyze the asymptotic behavior of 
a large class of total measures fn,ki our scheme of proof for the HGM system extends 
directly to (in this sense) arbitrary next neighbor interactions. This completes our 
presentation of the uniform case. 

3. Nonuniform case 

We now generalize the ideas of Section ^ to the nonuniform case. Our approach 
is formal. We assume that for a sufficiently large class of pair interactions V and 
external potentials U the corresponding free energy density functional exists and 
is convex. In fact, Chayes and Chayes || showed, that for our domain of interest 
of D = 1 particle systems with stable (see, e.g., |l2)) hard core next neighbor 
interactions, the set of external potentials for which a convex free energy density 
functional exists is nonempty. Thus, our calculus solves indeed a well-posed inverse 
problem for (in this sense) arbitrary next neighbor interactions. 

We construct our overcomplete free energy functional through a reduction of the 
original problem to a readily solvable multicomponent inverse problem. The blocks 
of i interacting particles [cf. for the uniform case the instances (2.4) respective 



( 2.22| )] are controlled by an effective local activity of component i of the mixture. 
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Hence, since the blocks interact by definition via hard core potentials, the interac- 
tion between components i and j of the mixture is represented by an effective pure 
hard core Boltzmann factor. Thus, we reduce the inverse problem for arbitrary 
next neighbor interactions to a type of polydisperse hard rod mixture problem. 
Consequently, we first extend the technique of proof of Vanderlick, Davis, and Per- 
cus jl3| for an additive mixture of pure hard cores in an external field to our case 
of polydispersity. Results for the local block densities are given next. Finally, the 
isomorphism between the original and our overcomplete description is established 
via a trivial identity of functional analysis. 

The presentation of the nonuniform case closes with the two examples already 
treated at the uniform level. The especially suggestive formal structure of these 
results provides evidence for our conjecture on a direct nonuniform isothermal- 
isobaric representation. 

3.1. Overcomplete free energy functional format. We start with the stan- 
dard description of nonuniform systems in direct form. So consider first the Gibbs 



measure (see, e.g., |l2J) of (1.1) at chemical potential /i (implicitly including the 
momentum contributions). Denote the local activity by z(x) ~ exp{/3[/i — U(x)]}. 
Then the total Gibbs measure becomes || [in Dirac notation ab = ab(y, x) — 
J R a(y,z)b(z,x) dz] 

E = l + Iz(I-ez)- 1 I, (3.1) 

where I is the identity, / represents the constant I{x) = 1, z refers to the matrix 
with elements z(y, x) — z(x)S(y — x), and e is given by ( |l.2| ). If huL^-too z(x) = 
sufficiently rapidly, then S exists. Hence, on following Percus J9|, |l~T|l , the local 
particle density g can be written as 

^ = aT^TT = ^S+(*Ms)S-(a:) (3.2) 

with 

S" = (I- ez)- 1 !, S+ = 1(1- ze) _1 , (3.3) 
and thus we have by (|l.2| ) as well as lim^^-too z(x) = the boundary conditions 

lim E~(x) = lim E + (x) = 1, (3.4a) 

X^ — OQ X — > + OC 

lim E~(x)= lim E + (x) = H, (3.4b) 

x— »-+oo X — oo 

lim q(x) = 0. (3.4c) 

x — >±oc 

The strategy is to determine z as a functional of g for (in the above sense) arbitrary 
V. The main stage to settle this inverse problem is the concept of overcomplete 
description, formulated next. 

Consider the Gibbs measure for a nonuniform multicomponent system on the 
line R |l^, O, where the local activity of component i is given by [henceforth a 
bar signifies a quantity within the overcomplete description] 

Wi(x, I) = z{fz)*-\x + \, x - |). (3.5) 

Here, x is the position of the block of i interacting particles, I + a its diameter 
with I being the distance between the boundary particles, and / is brought in via 
(O). Hence definition (p\5|), as the total isothermal measure of a i particle system 
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in a box (x — ^,x + j) with fixed boundary particles and Boltzmann factor /, 
controls the local aggregates. So the block description [cf. for the uniform case the 
examples (2.4) respective ( [2.22| )] is completed by the pure hard core Boltzmann 
factor between components i and j, given within the present parameterization by 

Eij(x,x' ,1,1') = hi + i> (x — x). (3.6) 

— +a 

Then, on assuming that the indices of the blocks are determined statistically Jl3|, |llj , 
the total Gibbs measure for such a mixture can be written as (cf. Jl3|, |TT| ]) 

Z = l + Fw{I-ew)- 1 I, (3.7) 

where e is the matrix with elements eij(x,x',l,l'), w represents Wij[x,x' ,1,1') = 
Wi(x, l)8(x — x')S(l — l')8ij, Ii(x) = 1 is the constant vector, and we used the Dirac- 
Einstein notation ab = aikbkj(y, x) = X^fcli Jr a ik(y, z)bkj{z, x) dz. A term by term 
analysis of the sums ([O]) and (3.7) shows 5 = H. Hence we conclude 

S = l + FiT)(I-ew)- 1 7 r , (3.8) 

which is the desired overcomplete description in direct form. Similarly, upon using 
Ei^it/ 2 )' 1 = (I — the equivalence of (3.7) with the Brannock-Percus 

format H can be formally established. Finally the local block densities are found 
to be 

Qi(x,l) = = l;E + (x,l)wi(x,l)E^(x,l) (3.9) 

olnwi{x, l) a 

with 

S- = {I-ew)- 1 !, H+ =F{l-we)- 1 . (3.10) 



Therefore by (3.6) and again \\m x ^± 00 z(x) = we have for I G R 

lim S~(a; ) /)= lim E + (x,l) = l, (3.11a) 

x — * — oo x — > + oo 

lim S _ (a;,/)= lim E + (x,l)=E, (3.11b) 

x — >-+oo x— > — oo 

lim Qi(x,l)=Q. (3.11c) 

x — >±oo 



Now, what makes the overcomplete description ( |3.8| ) more handy than the orig- 
inal formulation (|3.l| ) is, that as for an additive hard rod mixture e is of rank 1 on 
index space (cf. uW), and hence can be immediately written as 

ejj = hi±g5i> +a = Si+qhjt+q. (3.12) 

2 2 2 2 

Correspondingly, the technique of proof of |l3| for the free energy functional of 
additive mixtures of hard rods goes indeed through, with only minor modifications, 
for the case of our overcomplete free energy functional for arbitrary next neighbor 
interactions. 

The details are as follows. Introduce the one-dimensional auxiliary functions (cf. 



13|) 



i=l 



A (x) = \. / <5/'+a ( x ~ x')wi(x' , l')E (x' , I') d(x' , I'), (3.13a) 



2 



i=l 



^ + (x)=YJ Z + (x',l')w i (x',l')6 l , +a (x'-x)d(x',l'), (3.13b) 



2 
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as well as 

E-(x) = l+ I A-(y)dy, E+(x) = l + 



A + (y)dy. 



Then we find by fl3.12| ) and ( |3.10| ) the reduction (cf. jD| §) 



and via (3.9) that 



S ± (i J Z) = H ± (a:±i±2; 



A~(x) q (x) A + (x) g + (x) 



E E+(xY E E~(x) 
with the densities at the boundaries of the blocks of i interacting particles 



gf(x) = f Qi (x± l -^,l) 
Jm 



dl 



and 



Hence, on using the boundary conditions (3.11), we infer that (cf. jL3l) 



E+(x) dE-(x) 



Q (x), 



E~(x) dE+(x) 



E dx E dx 

If we substract these two equations, then we obtain 

1 d[E-(x)E+(x) 



g+(x). 



dx 



Q (x) - g + (x), 



which integrates upon using the boundary conditions ( |3.11 ) via ( |3.15 ) to 

g2(xg+(x) =l+ r [g -( y) _ e +( y)]dv = u ( x) . 



Thus we deduce from (|3.1l|) , fl3.15|) , ( |3.19| ), and ( |3.2l| ) that (cf. fl3|, §) 



In E (x) 



g (y) , , -0*0 

—r-r- dy, In - 



w(y) 



dy, 



and so we conclude by (3.11) via (3.15) that 



In 3 = 



g (y) 



dy = 



g + (y) 



dy, 



or more symmetrically 



InS 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



which generalizes the results (3.10) respective (4.30) of for dimerizing hard rods 
respective sticky cores to systems with arbitrary V. 
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Next, let us analyze the one-dimensional marginals gf of gi. According to (3.17) 
and Q, we have by ( |37To| ) and Q 



g-(x + § ) = ^H+(x, 0) • J z{fzf- x {x, x-l)E-(x- |, I) dl, 
Hence, on using the reparameterizations 



E + (x + ^l)z(fz) l - 1 {x + l,x)dl 



E+(y,0) (zfy- 1 z(y,x)dy, 



(3.25a) 
(3.25b) 

(3.26a) 
(3.26b) 

(3.27a) 
(3.27b) 



as well as ( 3.15 ), we infer the block densities equations 

q7(x) = ~5+(x).z(x-%). I (fzy- 1 (x-ly)E-(y-1)dy, 

Qt{x) = i [ E+iy + Dizff-^y^ + ^dyzix+^-E-ix). 
^ Jr 

In particular, we find for i — 1 the profile equation (cf. Q) 
Ql (x) = gf(xT%) = |s+(x + §)z(a;)S-(x- §) = «(ar)^(x + f ,ac - f ), (3.28) 

where the two-dimensional weight function ip (cf. j^, |ll|) evaluates by (3.21) and 
(p|) for y > x to (cf. {§) 



^(y^) = w (y) ex P 



y g~(*) 

x U{Z) 



dz 



lo(x) exp 



g+(z) 



ip(y,x) = yjuj(y)u){x) 



exp 



u(z) 
y g-(z) + g+(z) 



dz 



oj(z) 



dz 



(3.29) 



(3.30) 



Thus, upon substituting ( |3.28| ) into (3.27), and since gi is the density of the 
monomers, result (3.27) characterizes blocks as conformations of monomers cor- 
related by the attractive part of the interaction. This proposition, as the foregoing 
hard rod mixture form density functional (3.24), indicates a general pattern of 
transition from D = 1 to D > 1 (cf. Section |j)! 

Finally, to close the system of equations for the local particle density, we derive 
a subsidiary condition for the marginals g i in terms of g. Let us first rewrite ( |3.2[ ) 
in the form 



g(x) 



This yields 
1 



JlnS 

Slnz(x) 



E 



Slnz, ShiWi(y,l) 



d(y,l) 



5 In Wi(y, I) S\nz(x) 



2 ^MZXi^S/'O] 5lnz(x) 



d(y,l). (3.31) 



g(x) = - / Z+(y,l)(I-zf)- 1 (y+l,x)z(x)(I- fz)- 1 (x,y~l)Z-(y,l)d(y,l) 

(3.32) 
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via (|3.8|), (p.9[), and (p.5| ). Accordi ngly, we have through a reparameterization of 
the integral and by (|3.2|) as well as (3.15) that B 



(3.33a) 
(3.33b) 



S-(x-f)= / (I-jV)(z,y)H-(y)^, 



S+(z + §) 



E+(y) (I-zf)(y,x)dy, 



or that 



5~(x)= (I-fz)- 1 (x,y)E-(y-%)dy, 
E+(x) = [ Z+(y+1)(I-zf)- 1 (y,x)dy. 



(3.34a) 
(3.34b) 



Hence, on combining (3.34) with the block densities equations (3.27), we obtain 

1 



g-(x + f ) = =5+(.x + §)z(x)3-(x) 



1 , 



q + (x-1) = ^E + (x)z(x)E-(x-%), 

and so we conclude via Q, ( ^35| ), and ( jDs| ) that (cf. §) 

q+(x-$)q-(x + $) 



q(x) 



(3.35a) 
(3.35b) 

(3.36) 



which is the desire d clo sure relation. On starting from (|3.3l| ), the very simple 
formal structure of ( 3.36| ) is ce rtain ly due to the D = 1 next neighbor nature of our 
problem. Therefore, the form (3.31) has a higher dimensional analogue [as the result 
( 3.27 ) respective ( 3.24 )], whereas ( 3.36 ) has not. Proposition ( 3.36 ) completes the 
deduction of the nonuniform case. 

In passing, we point out that the main results of our analysis are completely 
characterized by the pressure-like quantities, termed effective pressures, p , defined 
as [||B1 

Q±(X) 



Pp ± (x) 



Lj{x) 



with the g ± brought in via ( p,18|) and lo by ( p.2l|) . We also introduce 

rv 



(3tt (y,x) 



ftp (z) dz. 



(3.37) 



(3.38) 



Upon using definition (3.37) and the corresponding results in the form derived so 
far, one can show that our overcomplete format free energy F for D — 1 particle 
systems with next neighbor interactions is given by 



3F\U] = min 
ee£>„ 



g(x) \fiU(x) + In z(x)\ dx — In E 



(3.39) 



with D n = {i £ L 1 (M) : t(x) > A Jmt(x) dx — n}, n fixed, the profile equation 



P[p-U{x)] = \n 6 i{x) --ln[i; + (x)i>-(x)}, 



(3.40) 
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where ^(x) = oj(x T § ) exp[- ( S 7 r ± (a; + § ,x - §)], V" = V> + by fl3.29| ), and the 
total Gibbs measure of hard rod mixture form (cf. p~4|) 



In 3= ^ / [f3p-(y)+pp + (y)]dy. 



(3.41) 



One further finds that the block densities equations transcribe to 

e - ( s ) = Qi (x) ■ J (f-^y (x-%,y+%)exp[-pn-(x-a,y)]dy t (3.42a) 

q+(x) = ^exp[-/37r+(y, x + a)] (|^/) ^(y - § ,ar + §) ■ e +(a:) (3.42b) 

with the subsidiary condition ( 3.36| ). These results extend — indirectly — the well- 
known (see, e.g., fl2| ) isothermal-isobaric technique from the uniform to the nonuni- 
form case. So we are led to conjecture that a direct nonuniform isothermal-isobaric 
approach should be possible as well. 

3.2. Examples. We apply the theory to the adhesive respective HGM system. 
This allows us to make contact with the associated results for the uniform case, 
Section ||, and to provide detailed evidence for our conjecture on the nonuniform 
isothermal-isobaric format. 

To present the examples most adequately, we give first closed expressions for the 
sums . Again, these simplifications are only possible for the D = 1 next neighbor 
case. Namely, if we start with ( 3.35| ), and use ( 3.33 ) as well as ( 3.28 ), we end up 
with 



g(s-f) 

Q~(x) 

g(s+f) 

Q+( X ) 



1+ / f(y + a- x)(3p + (y)exp[-f3Tr + (y,x)}dy, 



= 1+ / f{x + a-y)(3p (y) exp[-/?7r (x,y))dy 



(3.43a) 
(3.43b) 



via the subsidiary condition ( p. 36 ), (3.19), and through definition (3.37). This 
result is equivalent to the one-component case of formula (3.12) of §]. 



3.2.1. Adhesive interaction. Due to the contact nature of ( |2.3| ), we find by ( |3.43] ) 

l + XPp + (x), (3.44a) 



e(ar-f) 



Q (x) 

gQ + l) 

Q+(X) 



= 1 + X(3p (x), 



(3.44b) 



which are the generalizations of ( [2.17 ) to the nonuniform case. Consequently, we 
have that 



u)(x) = 1 



q(v) dy = 1 - t(x), 



and by (3.36) also the case i = 1 of (2.15) at the nonuniform level, 

q(x) 



0i 0*0 



[1 + X(3p- (x -§)][! + \(3p+(x+§)} ■ 



(3.45) 



(3.46) 
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Moreover, rewriting (3.44) via (3.37) and (3.45) as 

= pp-(x)[l+\0p + (x)], 



g(s-f) 



gives 



1 - t(x) 
Q(x + §) 
l-r(af) " 

t'{x) 



Pp+(x)[l + \pp-(x)} 



2A 



1 ± 



1 - t(x) 



4Acr(a;) 
1 - t(x) 



Xt'(x) 



1 - t(x) 



(3.47a) 
(3.47b) 

(3.48) 



with a(x) = \[q{x + f ) + q(x — §)]. Hence, on substituting ( |3.48 ) into (3.41), we 
arrive at formula (5.17) respective (6.16) of j^]. In the uniform limit, we recover 
( pTl8| ) respective ( pTTg| ). 

Finally, if we specialize (3.42) to ( |2.3| ), then we obtain via ( |3~45| ) that 

Qi(x-ja) 



gi(x)^g^(x)y- 1 l[- [ 



3 = 1 
t-1 



-(x) = et(x) \ ll Hj 



3 = 1 



- t(x - ja) ' 

gt(x + j a) 

- r(x+ja) ' 



(3.49a) 
(3.49b) 



which correspond immediately to ( p. 15 ) for i > 1 through (3.4£) and (3.4' 



3.2.2. HGM system. For an interaction characterized by (2.21), the block densities 
equations (3.42) are no longer analytically tractable. Thus, upon substituting (2.21) 
into (3.43), we are left with 



Q(x ~ f ) 
g-(x) 

e(z + f) 

g+(x) 

and so we have by ( |3.36| ) that 

g(x) 



= 1 + A{1 - exp[-/37r+(.T + d,x)}}, 



1 + A{1 — exp[— (3tt (x,x — d)]}, 



(3.50a) 
(3.50b) 



(l + \{l-exp[-l3n-(x-%,x-%-<I)]}) 
gi(x) 

x (l + A{l-exp[-/57r+(a; + f + d, x + §)]}). (3.51) 

These results generalize ( 2.27| ) respective (2.26) for i = 1 to the nonuniform case. 

Furthermore, on substituting (3.43) into (3.21), we find through ( ^2l| ), ( ^29| ), 
as well as definitions (3.37) respective (3.35) 



u>(x) = 1 — t(.t) — A 



x+d 



g + (y) 



(ip (z) exp[— /3n (y, z)] dz > dy, (3.52) 



y-d 



yielding via definition ( 3.38| ), ( 3.50| ), ( [3.2S| ), and definition ( |3.37[ ) 

rx+d e(» + f) 



Uj(x) = 1 — t{x) + 



r:.-, V {.h (//.;/- </)]( 1 + -) -1 



-dy 



x+d 



f3p + (y)uj(x) exp[-/37r+(y, x)} dy, 



(3.53) 
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and so we conclude on using again proposition (3.50) that the effective pressure p 
is determined by [cf. for the uniform case ( [2.30 )] 



gQc-f) g(y+f) , 

i-^-WWJ, 7 i> g ' ( } 

exp[/37r-(y,y-d)] 1 + - 



A 



and likewise p + by 



exp[/37r+(y + d, y)]l 1 + - 

which is equivalent to formula (6.15) of M. This completes our discussion of the 
nonuniform HGM system within the p ± representation. 

In summary, the examples demonstrate in detail that a suggestive formal corre- 
spondence between the usual uniform isothermal-isobaric method and the nonuni- 
form effective pressures format exists. This supports the idea of an isothermal- 
isobaric representation for arbitrary nonuniform D — 1 particle systems. 

4. Conclusion 

We presented a realization of the overcompleteness strategy for D = 1 systems 
with next neighbor interactions, where the supplementary variables are the local 
block densities gi of i associated particles. The corresponding free energy den- 
sity functional is of the simple additive hard rod mixture type. The isomorphism 
between our overcomplete format and the standard description was established 
through a trivial identity of functional analysis. The first results on this class of 
density functionals have been derived by Kierlik and Rosinberg M for the special 
case of dimerizing hard rods. So we lifted their approach to particle systems with 
arbitrary next neighbor interactions. Moreover, we showed that all our main results 
are completely characterized by the effective pressures p^ . This proposition in con- 
junction with the detailed results for the adhesive interaction and the HGM system 
support the idea of a nonuniform isothermal-isobaric representation (analoguous 
to the well-known uniform isothermal-isobaric method) for arbitrary nonuniform 
D = 1 systems. 

To introduce our concept of overcomplete description most concisely, we started 
with the uniform case. We treated first the two examples of adhesive systems and 
HGM systems. These results formed also the basis for our conjecture on a nonuni- 
form isothermal-isobaric representation. We deduced our overcomplete description 
for arbitrary next neighbor interactions by means of an application of large devia- 
tion techniques, accompanied by the Hardy-Ramanujan formula on the number of 
partitions of a positive integer. 

What we can learn from the D = 1 case is the logic of the problem for D > 1 
systems with finite range interactions. In abstracting from the block densities 
equations ( 3.42) ) , (i) we decompose a given D > 1 system into blocks of i monomers 



correlated by the attractive part of the interaction, (ii) Hence, on interpreting the 
blocks as components of a mixture, where the interaction between the components 
is given by construction by a hard core potential, we arrive at the extension of 



the hard rod mixture form density functional (3.41). (iii) Finally, we establish the 



isomorphism between our overcomplete description and the standard formalism 
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by an application of an identity of functional analysis as in (3.31). A familiar 
strategy to organize such a first-principles approach is to reduce the geometry of 
the phase space by going over to the corresponding discrete systems (see, e.g., ]l2|]). 
In particular, the prototypical simply connected infinite lattices, along the Cayley 
tree (for an overview see, e.g., |l(|), can be used to demonstrate the effectivity of our 
method. Progress in this direction will be the subject of a future communication. 

Complementary to the reduction of the geometry of the phase space, another 
possibility when trying to apply mathematical techniques is to restrict the class 
of interactions. For the prominent example of systems with strong, short range 
interactions, free energy models can be constructed along these lines. We begin 
with the physical observation, that in such, e.g., D = 3 systems the blocks take 
on predominantly elementary sphere-like or cube-like conformations. Thus, the 
available analytic expressions for the hard sphere respective hard cube mixture free 
energy functional can be used to model the hard core interaction between the blocks 
[cf. (ii)]. Modelling of the internal block free energy [cf. (i)] and the incorporation 
of the subsidiary condition [cf. (hi)] requires a more incisive approach, but is still 
amenable due to the strong, short range character of the interaction. Also the 
implementation of this program will be reported in the future. 

References 

[1] G. E. Andrews, Number theory, Dover, 1994. 

[2] G. R. Brannock and J. K. Percus, Wertheim cluster development of free energy functionals 
for general nearest-neighbor interactions in D = 1, The Journal of Chemical Physics 105 
(1996), 614. 

[3] J. T. Chayes and L. Chayes, On the validity of the inverse conjecture in classical density 

functional theory, Journal of Statistical Physics 36 (1984), 471. 
[4] R,. S. Ellis, Entropy, large deviations, and statistical mechanics, Springer, 1985. 
[5] K. F. Herzfeld and M. Goeppert-Mayer, On the states of aggregation, The Journal of Chemical 

Physics 2 (1934), 38. 

[6] P. Kasperkovitz and C. Tutschka, Cluster analysis of a one- dimensional continuous system 

with nearest-neighbor interaction, Ukrainian Journal of Physics 45 (2000), 488. 
[7] E. Kicrlik and M. L. Rosinberg, The classical fluid of associating hard rods in an external 

field, Journal of Statistical Physics 68 (1992), 1037. 
[8] J. K. Percus, Eguilibrium state of a classical fluid of hard rods in an external field, Journal 

of Statistical Physics 15 (1976), 505. 
[9] , One-dimensional classical fluid with nearest-neighbor interaction in arbitrary exter- 
nal field, Journal of Statistical Physics 28 (1982), 67. 
[10] , Density functional theory in the classical domain, Recent Developments and Appli- 
cations of Modern Density Functional Theory (J. M. Seminario, ed.), Elsevier, 1996, p. 151. 

[11] , Nonuniform classical fluid mixture in one- dimensional space with next neighbor 

interactions, Journal of Statistical Physics 89 (1997), 249. 
[12] D. Ruelle, Statistical mechanics: rigorous results, Benjamin, 1969. 

[13] T. K. Vanderlick, H. T. Davis, and J. K. Percus, The statistical mechanics of inhomogeneous 
hard rod mixtures, The Journal of Chemical Physics 91 (1989), 7136. 

[14] M. Q. Zhang, The nonuniform hard-rod fluid revisited, Journal of Statistical Physics 63 
(1991), 1191. 

Departamento de Matematicas, Universidad Carlos III de Madrid, Avenida de la 
Universidad 30, 28911 Leganes, Spain 



Departamento de Matematicas, Universidad Carlos III de Madrid, Avenida de la 
Universidad 30, 28911 Leganes, Spain 



